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A method for generating a mesoscopic superposition state of the collective spin variable of a gas of 
atoms is proposed. The state consists of a superposition of the atomic spins pointing in two slightly 
different directions. It is obtained by using off resonant light to carry out Quantum Non Demolition 
Measurements of the spins. The relevant experimental conditions, which require very dense atomic 
samples, can be realized with presently available techniques. Long-lived atomic superposition states 
may become useful as an off-line resource for quantum computing with otherwise linear operations. 
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Generation of superpositions of macroscopically differ- 
ent quantum states have been attracting fundamental in- 
terest for a long time Q since they illustrate the super- 
position principle in an extreme setting. Recently, such 
states have also been proposed as an off-line entangle- 
ment resource for a quantum computer, with which the 
computer can work using only linear operations Q| . Non- 
deterministically produced "cat" states of an ensemble of 
particles serve in this proposal as an alternative to opera- 
tions with maximally entangled quantum states in other 
quantum computing schemes. 

Superpositions involving a few microwave photons in a 
cavity have been demonstrated in Q . Quantum superpo- 
sitions have been recently also demonstrated using cur- 
rents in superconducting interference devices |^|. There 
are proposals for realising such states in different phys- 
ical systems, see for instance for those which involve 
Bose Einstein Condensates (BEC). Methods for gener- 
ating superpositions of two coherent states \a) + \ — a) 
of the electro-magnetic (e.m.) field have been suggested, 
but they require a perfect measurement of the photon 
number which is impractical at present. 

In this paper, we propose a method for realizing a su- 
perposition state of the collective spin variables of a gas of 
alkali atoms. Specifically the state that would be realised 
is a coherent superposition of the atomic spins pointing 
in two slightly different directions. Our method, in which 
the state of the atomic spins is generated via an interac- 
tion with light pulses, is another step within the growing 
field devoted to developing a quantum interface between 
light and atomic ensembles. Strong coupling between col- 
lective degrees of freedom of multi-atom ensembles and 
multi-photon pulses of light can be achieved in free space 
and has been recently used to demonstrate squeezing of 



the collective spin I nfill and entanglement between two 
atomic samples lfj.ll.2j. This should be contrasted with 
the strong coupling of individual photons to individual 
atoms, which requires high finesse cavities 0,0- 

The method we propose here for generating a coher- 
ent superposition of the spins advances these techniques 
beyond the continuous variables approximation. Indeed 
in our proposal a strong atom-photon coupling has to be 
realized with a small number of atoms. The relevant ex- 
perimental conditions, which require very dense atomic 
samples, seem to be realisable using presently available 
techniques. The transfer of the atomic ensemble spin 
state onto the state of light is also in principle possible 
0, ^| , and therefore superpositions of coherent states of 
light could also be realised using our method. 

The collective spin variables for the ground state of an 

atomic ensemble J = 2J. j , where j is the total an- 
gular momentum of the i'th atom in the ground state, 
obey the usual commutation relations [J x , J y ] = iJ z . For 
simplicity we assume that each atom is a spin 1/2 sys- 
tem. We start with the coherent spin state with maximal 
total angular momentum (i.e. all spins parallel) where- 
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(4j- + 1) where N a is the number 
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of atoms. We will suppose that the evolution of the col- 
lective spin is limited to unitary transfirmations without 
dissipation. Under these conditions the total length of 
the spin vector is conserved, i.e. it stays on the surface 
of the same Bloch sphere. Throughout the paper the spin 
vector stays polarized close to the +z direction (close to 
| tz) state), but not necessarily exactly in this direction. 

We can then define the rescaled variables xa — \ ~fs~ Jxi 
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which satisfy 



[xa,Pa] = i + 0(h A /N A ) 
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Thus xa and pa are analogous to the position and mo- 
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mentum operators, or to the quadrature operators of a 
single mode of the e.m. field, and fiA, the number of 
spins polarised in the +z direction, is analogous to the 
occupation number of a harmonic oscillator, or to the 
number of photons in a mode of the e.m. field. The new 
feature of the present model, as compared to the earlier 
work, is that it goes beyond the harmonic oscillator ap- 
proximation since the deviation of J z from N a /2 is taken 
into account to first order. 

The qualitative outline of the procedure for the "cat" 
state generation goes as follows, see figure ^ As shown 
in [E EH a component of the collective spin operator 
can be measured in a QND way. A measurement of J x 
to precision A J produces a squeezed state of the atomic 
spin with uncertainties Aa; 2 ^ = 2AJ 2 /N a = l/(2£ 2 ) and 
Ap^ = £ 2 /2 where £ denotes the degree of squeezing. 
Such a state has a "banana" shape localized around the 
pole of the Bloch sphere. The critical next step is to 
carry out a QND measurement of J z , i.e. of the number 
operator ua, with the same precision A J. Such a mea- 
surement can be visualized as intersecting the "banana" 
with a plane parallel to the equatorial plane ^ The dis- 
tance from the plane to the pole depends on the result of 
the QND measurement of ua- The "width" of the plane 
depends on the precision A J with which ua is measured. 
The resulting state lies at the intersection of the plane 
and the initial squeezed state. It is a superposition state. 
Of central importance is that only "weak" QND measure- 
ments of the atomic population, i.e. measurements of ua 
which do not need to have the single atom accuracy, are 
required to produce the superpostion state. 

Let us now describe in detail the QND measurements 
of the atomic spin. The starting point is the effec- 
tive interaction between the off resonant light and the 
atomic ground state spin 0, • Atomic levels and 
interactions with light are shown in figure la). Let 
us suppose that initially the photons are in the coher- 
ent state \ipp) = J dxp exp[— x%/2]\xp) (hereafter we 
omit for simplicity normalisation constants) and that 
the atoms are polarised in the z direction in the co- 
herent spin state \0) nA = / dxA exp [— x 2 A /2\ \xa) (the 
level scheme on top in figure la). For light propagat- 
ing along the x axis the effective interaction Hamilto- 
nian is H eff = nJ2i \ \a\a L + \ l x )i{ix \a R d R 

where cll, dp are the annihilation operators for left and 
right circular polarized light (the interaction is picto- 
rially shown in the middle level scheme in figure la). 
We introduce the Stokes operators S x — d^d^ — d^dp, 
S y = d^dp + d R dL, S z = — id^dp + iaLaL and as- 
sume that S z ~ N p is close to its maximal value (where 
Np is the number of photons in the beam). Then we 
can define the rescaled position-mo mentum like opera- 
tors S x = xpy/Np/2 , S y = pp/^/Np/2 which obey 
[xp,pp] — i. With this notation the unitary evolution 
can be written (up to an unimportant overall phase) as 

U = exp I —ia J x S x ) = exp (— iaxAXp) (3) 



where a = a^J N a N p /2. 

In the first step a spin squeezed state is obtained as fol- 
lows 0, 0, 0] After the interaction with light the state 
is, inthepp basis, U\ipp)\0) nA = J dxAdpp exp[— (aXA — 
pp) 2 /2] exp [— x A /2~\ \xa) \pp)- Upon measuring pp a 
particular outcome pp is obtained. The state of the 
atoms conditioned on the measurement outcome pp is 

4>a(jPp) = J dx A exp[-(ax A -pp) 2 /2] exp \_-x\j2\ \x A ) ■ 

This is a squeezed state. By a rotation one can center 
the state around pp = whereupon one has 

= dx A exp[-(a 2 + l)x A /2]\x A ) 

= J2^A)\n A ) (4) 

UA 

where c(tia) — if ua is odd, and equals c{ua) — 

( 2 , \ n A /2 > p 
2(^+1) J (n a/2)! When 71 A iS eVeil > With £, 2 = Ct 2 + l. 

The next step is generation of the cat state from the 
spin squeezed state by carrying out a QND measurement 
of J z . Towards this end we send off resonant right cir- 
cular polarised light in a coherent state with the mean 
number of photons N p along the z axis through the sam- 
ple (the bottom level scheme in figure la). The tem- 
poral evolution of the light and spins is then given by 

U = exp (—iaJ2i \ \®r®r)- We define a trans- 
lated creation operator by dp = y/N p + d' R and the 
translated quadrature operators ip = [a' R + d R )/y/2, 
pE = —i(a'p — o>'r)/V2 so that the unitary evolution 
takes the form U = exp [— iafiA (N p + y/2N p xp)] where 
we have dropped a term proportional to a R d' R . The 
term proportional to ^V p does not modify the state of 
the light, but it rotates the atomic spin. This term can 
be taken into account at later stages of the protocol and 
is omitted from now on. The part of the unitary trans- 
formation relevant for the QND measurement is therefore 
U = exp [— ipfiAXp] where (3 = a^j2N p . 

Initially light is in the coherent state ipp = 
J dxpexp [— £p/2] \xp) and the spin squeezed state of 
atoms is given by eq. Q ■ After the interaction of light 
and atoms, one measures the pp quadrature and finds an 
outcome pp. The state of the atoms conditional on this 
outcome is 

^a(pp) = ^2 c(n A ) exp [-((3n A - Pp) 2 /2] \n A ) ■ (5) 

TIA 

For suitable values of a,/3,pp, this is a cat 
state. In order to show this one should ex- 
press the state in the xa and pa basis. This 
can be done using the matrix elements (paI^a) — 
7:- 1 ^(2^n A \r 1/2 H nA (pA)exp[-p 2 A /2} where H jm <j v ) 
are the Hermite polynomials. A numerical calculation 
of iPa{pr) in the pa basis is given in figure 2 (solid line). 
To facilitate the analysis of the obtained state we de- 
rive approximate analytical expressions for iPa(j>r) in the 
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xa and pa basis. First note that using the Stirling ex- 
pansion, the coefficients c(tia) can be approximated by 

/ f 2_-,\™A/2 

c(nA) — ( J when ua is even. Inserting this 
expression into eq. JSJ) we find that 



iPa(pr) ^ Er. 



exp [-/3 2 (rM-//(pfl)) 2 /2] |n A )(6) 



where »(p R ) = f + 4ln ( ) ~ p -f . (7) 



The distribution of number states is thus a Gaussian with 
the mean ^{pr) and the variance 1//3. We now use the 
fact that the initial state is strongly squeezed in xa, see 
eq. ijljl. Thus we can use the approximation ha — Pa/^ 
in eq. 10 • This implies that the wave function in pa rep- 
resentation is centered around the two values ±y/2n(p R ) 
with the variance l/(y/2n(p R )f3). Hence the wave func- 
tion is 



\Hpr)) 



dpA C 



■ exp 



exp 



exp 



4/3 2 ^ 



cos 



IPA) (8) 
xa\/2/I | a; a) (9) 



where the second equality if obtained from the first by 
taking the Fourier transform, and the notations are sim- 
plified as fi — fi(jjji). Note that the Hermite polynomials 
with even index are even functions of pa which explains 
the + sign between the two terms in eq. (JSJ. As seen 
in figure 2, eq. © is a good approximation to the exact 
wave function (dotted line). 

To verify that the cat state has indeed been obtained, 
the measurements of spin components x a and pa should 
be carried out. The results of the pa measurement should 
be distributed according to the 2 broad peaks described 
by eq. JSJ. The condition for observation of two distinct 
peaks is that the width of each of them should be less 
than the distance between them. The xa measurement 
should exhibit an interference pattern due to the cosine 
in eq. J^J. This interference pattern is the proof that 
the two peaks found in the pa measurement are coher- 
ent, and therefore that a quantum superposition between 
atomic spins pointing in two different directions has been 
created. The condition for observation the interference 
is that the period of the cosine should be less than the 
width of the Gaussian in In fact both conditions 

lead to the same inequality, /i > This condition 

means that by the QND measurement of the population 
one has removed the components with small tia, i-e. the 
components near the origin pa = 0, from the squeezed 
state leaving a state with two peaks. 

There is a second condition which is that the average 
population jj, is limited by the degree of spin squeezing 
obtained in the first step of our protocol: fi < £ 2 . This 
condition has two origins. The first is that the probabil- 
ity of finding values of \i which exceed £ 2 is exponentially 
small. The second is that when this condition is violated 



it becomes impossible to measure the interference pat- 
tern in xa basis, and thereby to prove that the state so 
obtained is indeed a coherent superposition. 

Putting these two conditions together we obtain the 
condition (3£ 2 > 1 which relates the precision l/£ with 
which one can measure the xa quadrature to the preci- 
sion 1 /fi with which one can measure the number opera- 
tor fiA- Using the expressions for (3 and £ in terms of a, 
we can reexpress this condition as 



(10) 



Thus a superposition state can be produced if a high de- 
gree of squeezing can be realised with a small number of 
atoms. In what preceeds we have assumed that the quan- 
tum states are pure and that there is no noise present. 
Below we consider decoherence of the atomic state caused 
by the non-ideal QND measurements and describe a pos- 
sible experimental realization of the proposal. 

The off resonant probe light used for the QND mea- 
surements will cause (weak) real excitations of atoms 
leading to an incoherent spontaneous process |12|. The 
significance of this process can be assessed in terms of 
accessible experimental parameters. A first key param- 
eter is the resonant optical depth of the atomic sam- 
ple, kq = <rN a /A where a is the cross section. The 
forward scattering amplitude of light is proportional to 
^/(iA + 7/2) where 7 is the linewidth and A is the 
detuning with respect to resonance. Hence the optical 
depth at detuning A ^> 7 is ka = ^o7 2 /4A 2 and the 
phase shift (or polarisation rotation) of light traversing 
the sample is 8 a — K07/2A. The parameter a defined 
in eq. (PJ is the phase shift or polarisation rotation per 
atom, hence equal to a — 0A/N a . The degree of squeez- 
ing that can be obtained is £ — a J N a N p /2. Decoher- 
ence caused by the probe light can be described by the 
degree of optical depumping, rj — kaJV p /JV„ which can 
be interpreted as the probability that an atom makes an 
incoherent transition. All these parameters can be put 
together in the equation £ 2 = kqv/^ relating the degree 
of squeezing, the resonant optical depth and the degree 
of optical depumping. The degree of optical depumping 
must obey the condition 77 < l/£ 2 otherwise the coher- 
ence of the squeezed state or the superposition state will 
be destroyed. This condition follows from the fact that 
a fraction l/£ 2 of the spins are pointing in the —z di- 
rection and the other spins are pointing in the +z direc- 
tion. Hence on average Na/£, 2 atoms can decohere before 
the coherence of the superposition state, or of the par- 
ent squeezed state, is seriously affected. Thus the figure 
of merit for generation of squeezing and entanglement 
is the resonant optical depth of the sample which limits 
the degree of spin squeezing to £ 2 < \tJK$. Using these 
conditions and eq. (|10|l . one finds that the condition for 
realising a superposition state can be written as 



Kq 



> 47V 2/3 



fll) 



that is a large optical depth must be achieved with a 
small number of atoms. This means that sending light 
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through the sample along orthogonal directions for the 
two steps of the protocol is not optimal. Rather it seems 
more practical to confine the atoms to a cylinder of area 
A and length I, with light sent along the long axis of the 
cylinder, and to rotate the atomic spins between succes- 
sive measurements. 

The condition eq. can be achieved in, e.g., an al- 
kali metal BEC. Atomic density of the order of 10 15 cm -3 
for 4 • 10 5 atoms has been reported in for a cigar 
shaped condensate with the optical density along the 
long axis on the order of 10 4 . This is marginally close 
to the condition eq. ljTT|) . The degree of squeezing 
required for cat state generation under the above con- 
ditions is £ 2 50. An additional improvement in the 
efficiency of the spin measurement of a small number of 
atoms can be achieved by placing the atomic sample in 
a low-finesse symmetric standing wave cavity with the 
mirror transmission coefficients T = 1 — R. If the single 
pass absorption probability and the single pass phase ro- 
tation angle 9 are much smaller than T, then the phase 
shift/polarisation rotation of light exiting the cavity is 
9 C = 29 /T and the degree of optical depumping becomes 
rf = 2r//T. Thus the effect of the cavity can be summa- 
rized by introducing an effective cross section a c = 2a/T. 
The condition eq. (|TT1) for realisation of a mesoscopic su- 
perposition state therefore becomes 2kq/T > AN^ 3 - For 
a BEC sample with the density of 10 15 cm -3 placed in 
a low finess cavity with T = 5%, this condition is met 
with 10 3 atoms, which according to eq. (|10|l requires a 
feasible degree of spin squeezing f 2 w 10 for generation 
of the superposition state. Note that the life time of the 
coherent superposition state is close to the life time of 
the parent spin squeezed state, t c £ -2 where r c is 
the coherence time of the atomic ground state, which for 
a BEC is at least 100msec |l6j . 

Finally let us mention some additional experimental 
conditions. We have assumed above that the initial 



atomic spin state is perfectly polarized. In practice the 
degree of spin polarisation obtained by optical pumping 
m a vapor can reach 99% @, which should be sufficient 
to generate squeezing up to approximately £ 2 = 100. The 
degree of spin polarization in a BEC can be even higher 
because only certain magnetic states can be trapped. An- 
other challenge associated with the cylindrical configura- 
tion is that the atomic spins must be very precisely ro- 
tated between each QND measurement (using either opti- 
cal or radio frequency pulses). It is easy to check that the 
precision of this rotation must satisfy A9 < l/(£ 2 \fNl) 
which for the free space scenario with a BEC corresponds 
to the precision A9 ~ 3 x 10~ 5 . For the low finesse cavity 
scenario a lower precision of A9 — 1/300 is necessary. 

In summary we have proposed a method for generating 
a coherent superposition of states in which the collective 
spin of an atomic ensemble points in two slightly different 
directions. It appears to be realisable using present tech- 
nology by combining the methods used for QND mea- 
surements of collective atomic variables with a BEC in a 
bad cavity. Preliminary experimental work towards re- 
alisation of this proposal is under way. This mesoscopic 
superposition state uses the ground state atomic spins 
which means that it is relatively long lived. It can there- 
fore in principle be used as an offline non deterministic 
resource for quantum logic as described in Q. We will 
report on this in a future publication and in particular 
we will show that the methods described here can be 
used for purifying the entanglement between two distant 
atomic ensembles. 
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FIG. 1: Generation of the superposition spin state, a) Atomic 
levels and interaction with light. Both ground and excited 
states are spin 1/2 states. The top diagram shows the initial 
coherent spin state in z representation. The medium diagram 
shows the same state in x representation along with the vir- 
tual transitions used for generation of the spin squeezed state. 
The bottom diagram shows the spin squeezed state in z rep- 
resentation along with the interaction used for the cat state 
generation, b) The spin squeezed state generated by a QND 
measurement using the first pulse of light is represented by 
the "banana" shape. The following QND measurement of the 
number operator fiA corresponds to intersecting the "banana" 
with a plane of a finite width. Thus the final state has two 
components lying at positive and negative values of J y , i.e. it 
is a superposition state. 
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FIG. 2: The solid line shows the amplitude of the state 
4>a(pr), eq. @, in the p_\ basis for parameters £ 2 = 20, 
Pr/P = 7, P = 1/3. The dotted line shows the approximation 
eq. JHJ for the same values of the parameters. 



